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An  order  of  magnitude  model  is  presented  to  estimate  radiation  from 
detonation  of  conventional  explosives  in  an  attempt  to  predict 
frequency  bands  and  signal  levels  detected  by  other  investigators.  An 
earlier  model  describing  the  radiation  generated  by  explosions  has 
been  refined  to  include  the  contribution  of  the  heat  capacity  of  the 
detonation  products  and  the  temperature  dependence  of  the 
concentration  of  ionized  particles.  Relationships  are  established 
between  explosions  of  uncased  Composition  B,  the  radiation 
frequency  bands,  and  the  E-  and  B-field  amplitudes  as  a  function  of 
detection  distance.  The  model  considers  the  radiation  from  particles 
ionized  by  passage  of  the  shock  wave. 

A  comparison  of  the  calculated  radiation  with  thermal  and 
background  noise  estimates  shows  that  the  radiation  is  not 
detectable  above  the  background  radiation  even  for  large  explosives 
at  close  distances.  The  fact  that  radiation  has  been  observed 
indicates  either  that  the  assumptions  over-simplify  the  phenomenon 
or  that  the  primary  mechanism  of  radiation  production  has  been 
overlooked.  Ionized  particles  exist  10^  to  10^  times  longer  than  the 
time  to  accelerate  across  the  shock  wave,  during  which  the  particles 
could  produce  radiation  by  some  other  means.  Therefore,  it  is  likely 
that  the  model  greatly  underestimates  the  quantity  of  radiation 
produced. 
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1.  Introduction 


1.1  Background 

This  work  is  the  result  of  a  U.S.  Army  Research  Laboratory  (ARL)  technol¬ 
ogy  base  program  to  investigate  the  causes  of  electromagnetic  radiation  from 
the  detonation  of  conventional  explosives  in  support  of  future  combat  sys¬ 
tems.  The  overall  program  objective  was  to  provide  the  technology  to  de¬ 
velop  a  passive,  portable,  self-contained  unit  that  can  use  input  from  mul¬ 
tiple  sensors  to  detect,  locate,  and  identify  sources  of  weapons  fire  and 
explosives  detonation.  An  estimate  of  the  radio  frequency  (RF)  radiation 
from  the  detonation  of  conventional  explosives  is  presented  here,  which 
updates  an  earlier  ARL  model  [1]. 

Previous  reports  [1-3]  have  addressed  several  aspects  of  the  problem.  Ref¬ 
erence  [1]  provided  an  order-of-magnitude  model  based  on  the  thermal  ef¬ 
fects  of  explosive  materials  in  producing  ionized  particles  that  generate  the 
radiation.  Emphasis  was  on  blast  effects  that  had  propagated  far  from  the 
actual  detonation.  Reference  [2]  addressed  the  problem  of  RF  effects  on  a 
projectile  fuze  as  it  passed  directly  through  a  shaped  charge  explosion. 
Emphasis  was  on  effects  close  to  the  detonation  region,  and  information 
was  provided  about  the  heat  capacity  of  the  detonation  products  and  ef¬ 
fects  of  temperature  in  the  explosive  region  on  the  concentration  of  ionized 
particles.  Reference  [3]  was  a  literature  survey  of  previously  detected  RF 
energy  associated  with  battlefield  munitions  [4-12]. 

This  report  describes  refinements  of  the  model  of  Reference  [1]  that  include 
consideration  of  the  contribution  of  the  heat  capacity  of  the  detonation  prod¬ 
ucts  and  the  temperature  effect  on  the  concentration  of  ionized  particles. 
We  estimate  the  characteristics  of  the  RF  emission  from  a  detonation  of  a 
bare  explosive  charge  about  the  size  of  a  mortar  projectile. 

1.2  Summary  of  Literature  Survey 

The  most  prominent  results  of  earlier  work  are  stunmarized  in  table  1.  Pre¬ 
vious  investigators  detonated  explosives  including  lead  azide,  RDX  (rap¬ 
idly  detonating  explosive).  Composition  B,  Tetryl,  and  trinitrotolul  hexogen. 
Charge  sizes  from  fractions  of  a  gram  to  500  tons  (5  x  10^  grams)  were  re¬ 
ported.  Some  radiation  was  from  large  caliber  Navy  projectiles,  and  other 
radiation  was  from  bare  explosives  or  metal-encased  explosives.  Observed 
radiation  frequencies  were  reported  from  1  Hz  to  1  GHz.  Delay  times  from 
initiation  of  the  detonation  to  observation  of  the  radiation  ranged  from  80 
to  600  ps. 

Two  investigators  reported  E-field  values  at  distances  close  to  the  detona¬ 
tion.  Takakura  [5]  reports  4  x  10^  V/m  in  the  frequency  band  from  6  to 
90  MHz  at  distances  of  Im  from  0.1  to  0.4  g  of  lead  azide.  Using  the  inverse 
distance  dependence  for  the  radiated  E-field,  we  have  extrapolated  these  to 
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4  X 10"®  V/mat  10  km.  Curtis  [7]  reported  5  x  10^  V/mfrom  lOg  of  RDX  at 
a  distance  of  6.1  m  in  the  frequency  range  up  to  350  Hz.  This  extrapolates  to 
3  X  10"^  V/m  at  10  km. 

In  most  cases,  the  researchers  suggested  possible  physical  mechanisms  (see 
table  1)  but  presented  no  analytical  model.  Modeling  is  difficult  because 
there  may  be  multiple  causes.  Furthermore,  radiation  from  tube-launched 
weapons  such  as  small  arms,  mortar,  and  artillery  may  exhibit  radiation 
that  is  comparable  to  bomb  or  projectile  detonations  because  the  tube  walls 
conduct  charged  particles  and  concentrate  the  expanding  propellant  gases 
in  one  direction. 

Table  1.  Frequency  bands  observed  by  investigators  reviewed. 


Investigator 

Type  of 
explosive 
used 

Amount  of  Delay/  Frequency 

explosive  duration  of  range 

used  obser\^ed  signals 

Possible  cause 
suggested  by  authors 

Experimental  values  of  frequency  ranges 

Trinks 

Tube-launched 

artillery 

projectiles 

None  given 

1-100  kHz 

2  MHz-lGHz 

10  MHz-2GHz 

Muzzle  flash,  ionization 
of  gases  near  muzzle. 

Pulses  upon  impact  at 
target. 

Radiation  at  detonation 
from  "microsparks" 
caused  by  charge 
equalization  at 
detonation. 

Takakura 

Lead  azide 

0.1-0.4  g 

80-100  )is  delay 

6-90  MHz 

Acceleration  of  electrons 
ejected  by  ionization 
and  dipole  formation  at 
shock  front. 

Stuart 

Large  caliber 
guns 

— 

250  MHz-1  GHz 

None  given, 
experimental  results 
only. 

Curtis 

RDX 

lOG 

2-s  delay/19  s 
duration 

0.5-350  Hz 

None  given, 
experimental  results 
only. 

Gorshunov  et  a  I 

50/50 

trinitrotolul 

hexogen 

1000-5000  g 

3  Hz-20  MHz 

Electrical  charges 
generated 
asymmetrically 
from  scattered 
electrified  detonation 
products 

Cook 

Composition  B 

70-1100  g 

Below  10  kHz 

Gaseous  detonation 
products  form  a  plasma 
at  surface  of  gas  cloud 
from  ionization  by 
passing  through  earth's 
electric  field.  Gas  cloud 
discharges  on  contact 
with  ground. 
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Table  1  (cont'd).  Frequency  bands  observed  by  investigators  reviewed. 


Investigator 

Type  of 
explosi\'e 
used 

Amount  of 

explosive 

used 

Delay/  Frequency 

duration  of  range 

obser\^ed  signals 

Possible  cause 
suggested  by  authors 

Experimental 

yalues  of  frequency  ranges 

Wouters 

None  given 

1,300  g 

500  ton 
(=  4.5  X 

10=^  g) 

None  explicitly  — 

giyen 

8-ms  duration 
(1.3  kg) 

32-ms  duration 
(500  ton) 

Blast  temperature 
ionizes  detonation 
products  and  ambient 
air  and  produces  a 
plasma. 

van  Lint 

Bare  spheres 
to  metal- 
cased  bombs 

10-345,000  g  100-200-gs  50  MHz-1  GHz 

(bare  spheres  delay 

to  metal- 
encased  bombs) 

Separation  of  charge  at 
interface  of  explosion 
products  and  air  to  form 
a  vertical  dipole 
moment,  with 
asymmetry  induced  by 
reflection  of  shock  wave 
from  ground. 

Electric  sparks  from 
explosion  products 
interacting  with  casing 
fragments. 

Andersen  and 
Long 

Bare,  plaster- 
encased,  and 
seeded 
explosix’es 

Tetryl, 

Composition  B 

20-1 ,087  g 

300-600-|is  Less  than  600  kHz 

delay 

Detonation  ionizes 
detonation  products, 
which  transfer  charge  b\ 
friction  to  inert  casing 
particles  and  fragments. 
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2.  Theoretical  Model 


2.1  The  Detonation 


We  present  an  order-of-magnitude  estimate  of  the  electromagnetic  radia¬ 
tion  from  a  detonation  or  explosion.  We  use  the  term  ''explosion"  in  the 
sense  described  by  Wilfred  Baker  [13]  as  "a  process  by  which  a  pressure 
wave  of  finite  amplitude  is  generated  in  air  by  a  rapid  release  of  energy. 
Some  widely  different  types  of  energy  sources  can  produce  such  pressure 
waves  and  thus  be  classified  as  'explosives'  according  to  our  definition...." 

"Regardless  of  the  source  of  the  initial  finite  pressure  disturbance,  the  prop¬ 
erties  of  air  as  a  compressible  gas  will  cause  the  front  of  this  disturbance  to 
steepen  as  it  passes  through  the  air.. .until  it  exhibits  nearly  discontinuous 
increases  in  pressure,  density,  and  temperature.  The  resulting  shock  front 
moves  supersonically,  faster  than  sound  speed  in  the  air  ahead  of  it.  The  air 
particles  are  also  accelerated  by  the  passage  of  the  shock  front,  producing  a 
net  particle  velocity  in  the  direction  of  travel  of  the  front.  This  report  uses 
the  terms  "detonation"  and  "explosion"  interchangeably. 

The  model  is  used  to  estimate  the  frequency  bandwidth  of  the  radiation 
produced  and  the  electric  (E)  and  magnetic  (B)  field  amplitudes  versus  dis¬ 
tance  from  the  explosion.  Results  are  presented  at  10  km  from  the  explo¬ 
sion.  The  reader  can  obtain  results  at  other  distances  by  using  the  formulas 
that  are  provided.  The  model's  results  are  a  function  of  the  temperature 
within  an  explosive  region.  We  establish  this  temperature  dependence  by 
assuming  an  idealized  detonation  model,  is  shown  in  figure  1.  The  assump¬ 
tions  used  in  the  estimate  are  italicized. 


Some  of  these  assumptions  are 

•  The  iindetojiated  projectile  occupies  a  spherical  voliune  defined  by  radins  r^j. 

•  The  detonation  is  initiated  at  t=0  at  the  center  of  a  spherical  explosive  charge. 

•  After  initiation  of  the  explosion,  a  detonation  wave  propagates  outward 
through  the  explosive  charge  with  a  velocih/  of  10  ninilps  =  10^  in/s  [141. 

•  At  completion  of  the  detonation  at  t=^t^^  (3  ps  in  this  example),  the  detona¬ 
tion  wave  has  expanded  to  radius  r^y  and  all  the  iindetonated  explosive  has 
been  converted  to  detonation  products  with  the  release  of  thermal  energy 
into  the  volume  of  radius  r^^  [It.  lAV  neglect  the  portion  of  energi/  that  is 
converted  into  kinetic  energi/  of  bomb  fragments,  explosive  products,  and 
optical  and  acoustic  radiation. 

•  Ko  expansion  occurs  for  /</^^. 

The  temperature  of  the  region  with  r  <  Tq  at  f  =  /q  is  given  by 


T-l/o) 


C 


■+T 


V/Dd  Pro,/ 


lufb 


(1) 
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Figure  1.  Idealized 
detonation  model. 
(The  undetonated 
charge  occupies  the 
spherical  region  of 
radius ;  o-  After 
completion  of  the 
detonation,  a  shock 
wave  moves  out  with 
a  velocity  The 
thickness  of  the 
temperature 
discontinuity  across 
the  shock  wave  can  be 
calculated  by  the 
formulas  provided  in 
the  text.) 


Shock 

front 


Radiating  region 


'^amb 


in  which  =  300K  (ambient  temperature  at  sea  level), 


(2) 


and  C\  is  the  average  heat  capacity  of  the  material  inside  the  sphere  of  radius 

ro- 

For  t  >  tf,,  the  shock  loave  spreads  and  the  temperature  of  the  enlarging  spherical 
region  hound  by  r(t)  is  determined  by  the  heat  capacih/  of  the  additional  air  as  well 
as  the  detonation  products:  The  energi/,  temperature,  and  pressure  are  distributed 
uniformly  inside  the  spherical  region  behind  the  shock  wave. 

HV  assumed  that  the  specific  heat  capacih/  of  air  is  constant  and  equal  to  the  mean 
value  between  specific  heat  capacity  at  the  elevated  densihj  and  pressure  that  exists 
behind  the  shock  xoave  at  tg  and  specific  heat  capacih/  at  ambient  pressure  and 
temperature  ahead  of  the  shock  wave.  This  assumption  greatly  simplifies  our 
estimate  of  the  temperature  rise  inside  the  spherical  shock  region.  The  er¬ 
rors  as  large  as  a  factor  of  5  introduced  by  this  assumption  are  negligible 
compared  to  errors  of  several  orders  of  magnitude  introduced  by  the  tem¬ 
perature  variation  of  the  fraction  of  particles  ionized,  as  we  discuss  later. 

For  t  >  fg  and  r  <  r(t) 


T{t)  = 


^VAtiM  +  ^V/detproil 


+  T 


ninh  ' 


in 


which  C  V7,„v(f )  =  -  rf|)| 


Pnmb ^ \ s  • 


(3) , 

(4) 


Here,  is  the  density  of  air  at  ambient  temperature  and  pressure 
(1.23  kg/m  )  and  is  the  constant-volume  specific  heat  capacity  of  air  at 
ambient  temperature  and  pressure  (717  J /kg-K). 
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The  velocity  of  the  shock  wave  as  a  function  of  temperature  is  calculated  by 
Zel  dovich  and  Raizer  [15].  The  velocity  of  the  shock  wave  can  be  described 
by  their  results,  which  we  have  fitted  to  the  expression  (see  fig.  2): 

VJt)  =  5.5T"  m/s  for  (293  K  <  T  <  15,000  K) .  (5) 

Once  we  know  the  initial  velocity  of  the  shock  wave,  V^(to),  we  can  obtain 
the  radius  at  the  subsequent  times  by  using  the  following  equations: 

r(/,/  =  r(f,_,l  +  X  (/,-  f,  ,1  (6) 

in  which 

/,  =  /,_]  +  At  for  /  >  1 
with  At  =10’^  s. 

We  must  determine  the  value  of  C,^,,  in  equation  (3)  as  a  function  of 
temperature  and  pressure. 


2.2  Temperature  and  Pressure  Dependence  of  Heat  Capacity  of  Detonation  Products 

Composition  B  explosive  consists  of  63  percent  RDX,  36  percent  Trinitro¬ 
toluene  (dynamite)  (TNT),  and  1  percent  paraffin  filler.  The  model  assiime^^ 
that  the  detonation  tnke<  place  in  a  vacnnni,  i.e.,  in  the  absence  of  air.  Thus,  the 
program  does  not  include  further  combustion  of  the  detonation  products  in  air.  The 
combustion  product  gases  expand  outw^ard  adiabatically  after  the  comple¬ 
tion  of  detonation.  The  detonation  products  of  Composition  B  were  obtained 
from  the  Cheetah  computer  program  [16]  that  used  the  applicable  chemis¬ 
try.  The  detonation  products  obtained  from  the  Cheetah  model  in  concen¬ 
trations  greater  than  1  mole  of  detonation  product  per  kilogram  of  explo¬ 
sive  are  by  weight,  24.1  percent  H2O,  36.4  percent  N-,,  21.9  percent  CO2, 
5.2  percent  CO,  and  12.1  percent  solid  carbon. 


Figure  2.  Shock  wave 
velocity  as  a  function 
of  ten-iperature  behind 
the  shock  wave. 


The  specific  heat  capacity  of  the  detonation  product  gases  depends  on  tem¬ 
perature  and  pressure  [17],  [18].  Temperature  dependence  [17]  is  shown  in 
figure  3.  The  total  heat  capacity,  calculated  from  figure  3  and  the  relative 
concentrations,  varies  from  0.26  cal/gm-K  at  300  K  to  0.39  cal/gm-K  at  3500 
K.  The  pressure  dependence  is  given  in  [17].  The  !>pecific  heat  capacity  of  the 


Temperature  behind  shock  wave  (K) 
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Figure  3.  Temperature 
dependence  of  the 
specific  heat  capacity 
of  the  detonation 
products  of  Composi¬ 
tion  B.  (The  specific 
heat  capacity  of  each 
of  the  five  highest 
concentration 
detonation  products  of 
Composition  B  is 
shown.) 


Temperature  (K) 

detonation  products  varies  adiabatically  from  0.210  cal/gm-K  at  1  atm  to 
0.587  cal/gm-k  at  24  x  10* atm.  In  the  calculation,  we  use  the  mean  value  of  0.40 
over  these  temperature  and  pressure  ranges. 

The  pressure  and  temperature  variations  on  the  specific  heat  capacity  of  air 
vary  from  0.1719  cal/gm-K  at  1  atm  and  300  K  to  0.315  cal/gm-K  at  100  atm 
and  3000  K.  We  use  the  mean  value  over  this  temperature  and  pressure  range, 
0.243  cal/gm-K.  This  could  result  in  values  that  are  too  large  or  too  small  by 
30  percent  over  the  pressure  and  temperature  extremes. 

These  assumed  values  for  the  specific  heat  capacities  of  the  detonation  prod¬ 
ucts  and  ambient  air  lead  to  errors  that  are  negligible  compared  to  the  er¬ 
rors  arising  from  the  temperature  dependence  of  the  concentration  of  ion¬ 
ized  particles. 


2.3  Ionization  Fraction 

A  previous  report  [2]  showed  that  only  a  very  small  fraction  of  the  particles 
are  ionized.  We  therefore  need  a  method  that  accounts  for  the  ionization 
fraction  of  the  radiation.  The  Hilsenrath  and  Klein  calculation  [19]  of  the 
ionization  fraction  of  air  as  a  function  of  temperature  is  shown  in  figure  4. 
One  sees  that  the  ionization  fraction  is  very  sensitive  to  changes  in  tem¬ 
perature  below  15,000  K.  We  have  interpolated  the  lower  temperature 
ionization  fraction  results  from  Hilsenrath  and  Klein  in  figure  4.  live  inter¬ 
polation  formula  is 

Log^^IonizationFraction  =  Log^^fi  =  (-30.566)  -i-  (0.011267)  *  T  -  (l.3664x  lO"^)"^  (7) 

Table  2  shows  the  ionization  of  various  ion  species  in  air  at  3500  K  and  2500 
atmospheres  [19]  at  the  high  temperatures  and  pressures  that  could  occur 
during  an  explosion.  Since  air  is  composed  primarily  of  N2  and  O2,  and 
since  the  dominant  detonation  products  of  Comp  B  explosive  are  H2O,  N2, 
C,  O2,  and  CO2,  one  might  expect  ionization  fraction  values  of  the  same 
order  of  magnitude  as  air,  which  are  calculated  from  equation  7  and  plotted 
in  figure  4.  We  see,  however,  that  the  tabulated  results  from  table  2  are  quite 
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Figure  4.  Curve  fit  to 
ionization  fraction,  IF, 
of  air  as  a  function  of 
temperature.  (The 
logarithm  of  the 
ionization  fraction 
data  from  Reference  19 
is  plotted  versus 
temperature  of  the  air. 
The  line  obtained  from 
the  curve  fit  equation 
is  also  shown.) 


1500  2000  2500  3000  3500  4000  4500 


Temperature  inside  spherical  shock  region  (K) 


Table  2.  Ionization  of 
air  at  3500  K  and  2500 
atm. 


Species 

Ionization 

fraction 

Trend  with 
increasing  pressure 

NO* 

3.1  X 10-" 

decreasing 

O2* 

1.5  X  10-*' 

decreasing 

w 

2.3  X  10-’= 

decreasing 

0* 

2.5  X  10-’5 

decreasing 

CO* 

1.1  X  10-’^ 

decreasing 

E* 

1.0  X  10-« 

decreasing 

0* 

2.0  X  10^ 

increasing 

1.1  X 10’= 

increasing 

different  from  the  values  for  air  at  3500  K  in  figure  4.  Nevertheless,  we  assume 
that  we  can  use  equation  7  for  determining  the  ionization  fraction  of  the  detonation 
products.  We  further  assume  that  (a)  the  detonation  products  and  the  air  are  at 
most  singly  ionized,  and  that  (b)  the  electron  concentration,  e-,  is  equal  to  the  ion 
concentration. 


2.4  Radiation  From  Particle  Acceleration  Across  a  Shock  Wave 

Radiation  is  caused  by  the  acceleration  of  charged  particles  because  of  the 
passage  of  the  shock  wave.  Additional  assumptions  are  as  follow:  All  the 
acceleration  takes  place  over  the  thickness  of  the  shock  wave,  which  we  shall  shoiu  is 
on  the  order  of  10'^  m.  We  apply  the  equipartition  theorem  [201  to  provide  a  rela¬ 
tionship  for  the  mean  particle  speed  V.  Ideal  gas  behavior  applies,  and  the  air  before 
and  after  the  shock  wave  is  in  thermodynamic  equilibrium. 

Since  we  are  considering  only  the  radial  component  of  the  electron's  veloc¬ 
ity,  the  problem  is  one-dimensional,  with  the  result  that  the  speed  change, 
AV,  across  the  shock  (neglecting  direction  changes)  and  the  average  speed 
across  the  shock,  V,  as  shown  in  figure  5,  are  given  respectively  by 

=  ,  (8) 
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Figure  5.  Change  in 
mean  speed  of 
charged  particles 
caused  by  temperature 
change  across  the 
shock  wave  for  the 
one-dimensional  case 
considered  in  the 
analysis. 


in  which  k  is  the  Boltzmann  constant,  1.38  x  10"^  J /K. 

F(T) + v(r„,)  _  ) 

2  V  m  2  ' 


in  which  V(T)  is  the  average  radicd  velocity  for  r  >  Vq. 
The  shock  wave  Mach  number  is  defined  by 

y,m) 


M(0- 


V 


amb 


(9) 


(10) 


in  which  Vamb,  the  speed  of  sound  at  ambient  conditions,  =  340  m/s. 

The  mean  free  path  A,  which  is  the  average  distance  that  particles  travel  between 
collisions  [21],  is 


16  V 
5  -flltRTg^l, 


(11) 


in  which  v  is  the  kinematic  viscosity  of  air  at  ambient  conditions  (14  xl0~^  m^/s) 
and  R  is  the  ideal  gas  constant  of  air  (287  J/kg-K). 

The  shock  thickness  S(t)  is  given  by  Kogan  [21] 
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(5(0=^-|M(f)A=10-«M(0  (12) 

in  which  /is  the  ratio  of  the  constant-pressure  and  constant-volume  spe¬ 
cific  heat  capacities  of  air,  which  we  have  assumed  to  be  constant  over  the 
range  of  temperatures  and  pressures  behind  the  shock. 

The  volume  of  air  in  the  thin  shock  wave  is 


(13) 


With  the  assumption  of  constant-volume  and  constant-density  conditions,  the  pres¬ 
sure  rise  produced  is  the  bomb-yield  energy  divided  by  the  volume  of  the 
spherical  shock  region,  so  that  the  pressure  is 


(14) 


in  which  pg„i,  is  ambient  pressure,  1.01  x  10^  N/m^. 

During  these  conditions,  the  number  of  moles  of  air  in  the  region  in  the 
shock  wave  is  given  by  [22] 


n{t)  = 


p{t)Tjt) 

RT{t) 


(15) 


Now  consider  a  charged  particle  that  accelerates  with  an  acceleration  a (f) 
across  the  shock  front  for  a  short  time  At  and  continues  moving  in  the  same 
direction  with  the  new  velocity.  The  time  of  transit  of  the  charged  particle 
across  the  shock  wave  is 


and  the  acceleration  of  the  charged  particle  is  given  by 

At  S' 


(16) 

(17) 


We  assume  that  the  radiation  is  more  or  less  uniformly  distributed  over  a  band 
starting  at  0  Hz  and  extending  to  1/At  [23],  given  by 


^f^Af  (18)' 

and  that  from  equation  (17),  we  see  that 

^  ■  (19) 

The  time-average  total  power  per  singly  ionized  particle,  P,  integrated  over  all  di¬ 
rections,  is  given  by  [24] 


P  = 


eV 


_  _  Poefaf 
6tce(f^  ^  (>7t 


(20) 


in  which  e  is  the  electron  charge  (1.6  x  10’^^  C);  the  permeability  of  free 
space  (taken  as  equal  to  the  permeability  of  air),  is  1.26  pH/m;  amd  Sq,  the 
permittivity  of  free  space  (taken  as  equal  to  the  permittivity  of  air),  is 
8.85  pF/m. 
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The  total  power  radiated  for  n  moles  of  singly  ionized  particles  is 

P total  • 


(21) 


in  which  A  is  the  Avogadro's  number,  6.02  x  10^^  particles/ mole,  and/)  is  the 
fraction  of  particles  ionized  from  equation  (7). 

The  average  E-field  per  cycle  produced  in  the  far  field  per  particle  is  given 
by  [25] 


E 


ufYM 

ry  In 


(22) 


Using  equations  20,  21,  and  22,  we  obtain  the  following  formulas  for  the 
total  E-field  and  B-field: 


^total^  I  V 


(23) 


(24) 


Calculations  have  been  performed  for  the  case  of  an  explosion  having  a 
yield  that  approximates  the  energy  from  a  mortar  projectile  that  is  25  MJ. 


The  calculated  frequency  bands  shown  in  figure  6  include  the  observed  fre¬ 
quency  ranges  of  6  to  250  MHz  to  2  GHz,  as  shown  in  table  1.  The  electron 
frequency  band  is  1.4  THz,  which  is  too  high  to  be  observed  with  RF,  non- 
line-of-sight  equipment.  Figure  6  shows  that  the  upper  limits  of  the  fre¬ 
quency  bandwidths  for  N2  and  O2  ions  are  about  8  to  9  GHz. 

Figure  7  is  the  calculated  E-field  amplitude  at  a  range  of  10  km  versus  time 
from  initiation  of  detonation  for  the  O2  ions,  N2  ions,  and  electrons.  The  £- 
field  amplitude  depends  on  the  mass  of  the  particles  through  the  accelera¬ 
tion  (eqs  17  and  23),  the  velocity  change  and  average  velocities  (eqs  8  and 


Figure  6,  Calculated 
frequency  band  of 
radiation  from 
electron,  Oj  ,  and  N2 
ions  formed  when  the 
shock  wave  from  the 
detonation  passes 
through  ambient  air. 
(The  frequency  band 
in  volts /meter  is 
plotted  versus  time  in 
microseconds  from 
initiation  of  detonation 
at  10  km.) 
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Figure  7.  Calculated 
electric  field  ampli¬ 
tude  of  radiation  from 
electron,  O2  ,  and  A/ 2 
ions  formed  when  the 
shock  wave  from  the 
detonation  passes 
through  ambient  air. 
(The  E-field  amplitude 
in  volts/ meter  is 
plotted  versus  time  in 
microseconds  from 
initiation  of  the 
detonation  at  10  km.) 


9).  If  equations  (8)  and  (9)  are  substituted  into  equation  (17),  and  the  result¬ 
ing  acceleration  is  substituted  into  equation  (23),  the  resulting  expression 
shows  that  the  £-field  amplitude  is  inversely  proportional  to  the  masses  of 
the  particles.  Thus,  the  radiation  from  the  electrons  should  be  more  detect¬ 
able  than  the  radiation  from  the  heavier  0\  and 

If  we  arbitrarily  assume  that  the  detection  threshold  is  10“^  V /m,  then  from 
figure  7,  the  electron's  radiation  is  detectable  at  10  km  for  100  |xs,  but  the 
radiation  from  the  O2  and  ^2  ioris  is  not  detectable.  At  10  m  from  the 
explosion,  the  signal  from  the  ions  would  increase  to  2  x  10"^  V/m,  and 
would  remain  above  the  detection  threshold  for  about  3  ps. 


Time  from  initiation  of  detonation  (ps) 


12 


3.  Consequences  of  the  Model 


3.1  Effect  of  Explosive  Mass 

The  model  shows  that  the  mass  of  the  imdetonated  explosive  affects  both 
the  amplitude  of  the  radiation  and  the  delay  time  from  initiation  of  the  ex¬ 
plosion  to  completion  of  the  detonation. 

Equation  (23)  shows  that  the  E-field  amplitude  is  proportional  to  the  square 
root  of  the  number  of  moles  of  ionized  particles,  which  we  assume  is  propor¬ 
tional  to  the  mass  of  the  undetonated  explosive.  Thus,  the  model  leads  us  to 
expect  that  the  amplitude  ratio  should  be  proportional  to  the  square  root  of 
the  explosive  mass. 

The  delay  time  from  initiation  to  completion  of  the  detonation  is  approxi¬ 
mately  the  time  for  the  detonation  wave  to  propagate  at  constant  velocity, 
V^et,  from  the  center  of  the  spherical  explosive  (point  of  initiation)  to  the 
surface  of  the  explosive.  From  reference  [3],  the  ratio  of  the  delay  times  for 
two  spherical  explosives  of  uniform  density  but  different  masses  should  be 
proportional  to  the  ratio  of  the  cube  roots  of  their  masses. 

3.2  Radiation  From  Recombination  of  Ions  with  Electrons 

Heretofore,  we  have  discussed  the  radiation  caused  by  particles  ionized 
during  their  passage  through  the  shock  region;  we  have  not  considered  that 
the  ions  eventually  re-combine  with  electrons,  producing  further  radiation. 
Table  3  [26]  gives  the  ionization  energies  and  frequencies  of  the  detonation 
products  N2  and  CO2. 

The  frequencies  are  in  the  terahertz  region,  with  wavelengths  in  the  ultra¬ 
violet  spectrum,  and  therefore  are  not  useful  for  transmission  in  non- 
line-of-sight  scenarios. 

3.3  Detectability  of  RF  Signals  From  Detonation  of  Conventional  Explosives 

The  detectability  of  the  radiation  depends  on  the  thermal  noise,  background 
radiation  noise,  and  the  noise  inherent  in  the  sensor  and  data  acquisition 
system.  The  examination  of  particular  instrumentation  systems  and  back¬ 
ground  radiation  noise  is  beyond  the  scope  of  this  report;  however,  we  can 
make  some  general  observations  about  the  thermal  noise.  The  thermal  noise 
power  is  given  by  [27] 

^noise  ~  ^T'amb^f  (25) 


Detonation  Ionization  energy  Frequency  Wavelength 
product  (e-V)  (J)  (Hz)  (A) 


N2  15.5  2.48  X 10-“  3.73  x  10“  800 

CO2  14.4  2.30  X 10-“  4.0  x  10“  750 


Table  3.  Radiation 
from  recombination  of 
detonation  products. 
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in  which  is  the  absolute  temperature  ot  the  circuits  in  K,  here  taken  to 
be  300  K,  and  A/  is  the  frequency  bandwidth  of  the  signal  in  hertz.  The 
calculated  values  of  A/ and  E  from  figures  6  and  7  are  given  in  table  4,  with 
the  noise  power  calculated  from  equation  (25). 

The  power  Pg^,g  received  by  a  1-m^  antenna  far  from  the  source  can  be  esti¬ 
mated  from  the  plane  wave  relationship  [28],  via  E-field  values  obtained 
from  figure  7: 


|iceoEp.«A-pe#("2')  =  [0.00133E 
c  and  Co  were  defined  previously. 


P  „j,g{ioatts)  = 


2 

peak 


^eff 


(26) 


Table  4  shows  that  the  average  power  received  is  far  less  than  the  kTAf  noise 
as  close  as  1  m  to  the  detonation.  The  electron,  which  is  the  most  detectable 
particle,  is  11  dB  below  the  noise,  even  at  1  m,  and  is  91  dB  below  noise  at 
10  km.  Even  if  an  explosive  10,000  times  larger  than  the  mortar  (equivalent 
to  a  4000-lb  bomb)  is  used,  only  the  electron  radiation  is  detectable  above 
the  noise  at  close  distances. 


Using  E-field  amplitudes  and  frequency  bands  from  figures  6  and  7,  we  can 
calculate  the  signal  level  relative  to  the  noise  at  1  m  and  10  km  from  the 
detonation.  The  data  are  then  extrapolated  to  an  explosive  having  10,000 
times  more  mass. 

These  calculations  have  shown  that  the  wide-band  signals  can  be  expected 
from  radiation  caused  by  ionized  particles  accelerating  across  the  shock 
wave;  however,  the  signals  are  not  detectable  above  /cTA/ noise  for  normal 
sized  explosives  at  useful  detection  ranges.  Thus,  the  assumptions  have  over¬ 
simplified  the  phenomenon,  or  the  model  has  not  established  the  mecha¬ 
nism  responsible  for  the  signals  measured  by  previous  investigators  (see 
table  1  and  sect.  2). 


Table  4.  Detectability  of  estimated  radiation  from  shock  wave. 


Particle 

Calculated 

E-field 

(V/m) 

A/ 

(Hz) 

Noise  power 
(W) 

Calculated 
Average  power 
at  receiver 
(W) 

Detectability 

(dB) 

0.191  kg  explosive  (mortar) 

At  1  m 

Electron 

6  X  10-^ 

1.6  X  10'2 

6  X  10’ 

4.8  X  10-'" 

-11 

N2  Ion 

1  X  10-" 

9x10’ 

3.7  X  10-" 

1.3  X  10-'’ 

-84 

O2  Ion 

9  X  10-'’ 

8x10’ 

3.3  X  10-" 

1.07x10-” 

-85 

At  10  km 

Electron 

6  X  10-" 

1.6  X  10'2 

6  X  10’ 

4.8  X  10'" 

-91 

N2  Ion 

1  X  10-12 

9x10’ 

3.7  X  10-" 

1.3  X  10-22 

-164 

O2  Ion 

9  X  10-'" 

8x10’ 

3.3  X  10-" 

1.07x10-22 

-165 
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3.4  Discussion 


The  calculations  and  order-of-magnitude  model  address  the  radiation  pro¬ 
duced  by  explosions  that  occur  in  bomb  and  projectile  detonations.  The 
bare  explosive  (not  the  explosive  encased  by  the  bomb  or  projectile  shell)  is 
considered.  The  mechanism  of  radiation  analyzed  was  radiation  from  par¬ 
ticles  ionized  by  passage  of  the  shock  wave.  The  most  common  charged 
particles  from  the  atmosphere  are  the  O2  ions  and  Nj  ions  and  the  elec¬ 
trons.  The  calculations  that  were  presented  for  each  type  of  particle  are  the 
radiation  bandwidth,  time-average  E-field  and  B-field  amplitudes. 

The  values  of  bandwidth  measured  by  experimeters  (see  table  1)  ranged 
from  1  MHz  to  1  GHz.  The  values  from  the  model  extend  from  8  GHz  for 
the  ions  to  terahertz  for  the  electrons. 

The  heat  capacity  of  the  detonation  products  and  of  the  ambient  air  deter¬ 
mine  the  temperature  of  the  region  of  high  temperature,  pressure,  and  den¬ 
sity  behind  the  shock  wave.  The  concentration  of  ionized  particles  strongly 
depends  on  the  temperature  behind  the  spherical  shock  wave.  It  decreases 
very  steeply  as  the  temperature  drops  during  the  shock  wave  expansion. 
An  estimate  was  made  of  the  dependence  of  the  ion  concentration  on  tem¬ 
perature.  Relationships  were  established  between  explosive  mass  of  uncased 
Composition  B,  the  radiation  frequency  bands,  and  the  £-  and  B-field  am¬ 
plitudes  as  a  function  of  detection  distance.  At  a  range  of  10  km,  radiation 
from  free  electrons  caused  by  a  mortar-sized  blast  is  expected  in  the  fre¬ 
quency  band  up  to  1.6  THz  with  amplitudes  above  0.01  microvolt/m  for 
the  first  100  fis  from  detonation.  Oxygen  and  nitrogen  ions  produce  radia¬ 
tion  in  the  frequency  band  up  to  9  GHz  with  amplitudes  above  10  picovolt /m 
for  3  |is  after  detonation. 

A  comparison  of  the  calculated  radiation  with  thermal  noise  estimates  shows 
ikat  most  of  the  radiation  is  not  detectable  above  the  noise  even  for  large 
explosives  at  close  distances.  The  fact  that  signals  have  been  observed  sug¬ 
gest  either  that  the  assumptions  over-simplify  the  phenomenon  or  that  the 
radiation  detected  must  be  caused  by  other  mechanisms. 

The  times  for  the  ions  and  electrons  to  cross  the  shock  wave  are  the  recipro¬ 
cals  of  the  frequency  bands  from  figure  6: 10'^°  s  for  the  ions  and  lO'i^  s  for 
the  electrons.  The  particles  exist  in  the  charged  state  for  the  duration  of  the 
explosive  event,  which  from  figures  6  and  7,  is  nominally  100  ps  or  10^  s. 
This  is  10®  times  longer  than  the  shock  transit  time  for  the  ions  and  10®  times 
longer  for  the  electrons.  During  this  much  longer  time,  the  charged  par¬ 
ticles  could  be  producing  radiation  by  some  other  mechanism.  We  there¬ 
fore  believe  that  the  model  greatly  underestimates  the  quantity  of  radiation 
produced. 
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